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ABSTRACT 

Aims. We study the dependence of protoplanetary disk evolution on stellar mass using a large sample of young stellar objects in 
nearby young star-forming regions. 

Methods. We update the protoplanetary disk fractions presented in our recent work (paper I of this series) derived for 22 nearby 
(<500 pc) associations between 1 and 100 Myr. We use a subsample of 1428 spectroscopically confirmed members to study the 
impact of stellar mass on protoplanetary disk evolution. We divide this sample into two stellar mass bins (2 M 0 boundary) and two 
age bins (3 Myr boundary), and use infrared excesses over the photospheric emission to classify objects in three groups: protoplanetary 
disks, evolved disks, and diskless. The homogeneous analysis and bias corrections allow for a statistically significant inter-comparison 
of the obtained results. 

Results. We find robust statistical evidence of disk evolution dependence with stellar mass. Our results, combined with previous 
studies on disk evolution, confirm that protoplanetary disks evolve faster and/or earlier around high-mass (> 2 M G ) stars. We also find 
a roughly constant level of evolved disks throughout the whole age and stellar mass spectra. 

Conclusions. We conclude that protoplanetary disk evolution depends on stellar mass. Such a dependence could have important 
implications for gas giant planet formation and migration, and could contribute to explaining the apparent paucity of hot Jupiters 
around high-mass stars. 

Key words. Planetary systems: protoplanetary disks, planet-disk interactions - stars: formation - (stars:) planetary systems - stars: 
pre-main-sequence 


1. Introduction 

Studying the evolution of protoplanetary disks around young 
stellar objects (YSOs) is crucial to better comprehend how plan¬ 
ets form. Understanding the effects of fundamental parameters 
such as age, stellar mass, environment, or binarity provides in¬ 
sight into the planet formation processes, and could shed light 
onto important relations between planetary systems and their 
stellar hosts. Several studies on this topic (see e.g. Williams & 


Cieza||20TTj [Alexander et al. |2013| |Espaillat et al. 2014| ) have 


already drawn a comprehensive picture of protoplanetary disk 
evolution. 

In the first paper of this series ( |Ribas et al.||2014| hereafter 
R14) we presented a large sample of ~2 300 spectroscopically 
confirmed YSOs in 22 young (< 100 Myr) and nearby (< 500 pc) 


The evolution of protoplanetary disks and planets are closely 
connected as the disk properties define the initial conditions of 
planet formation, and the planet evolution subsequently affects 
the disk properties (e.g. Kley & Nelson 2012; Baruteau et al. 
|2013| >. The disk lifetime is a fundamental parameter to under¬ 
stand the interplay between a disk and its planets, since it deter¬ 
mines the time available for the planets to form and migrate (e.g. 
|Burkert & Ida |2007[ [Alexander & Armitage|[2009] ). A depen¬ 
dence of the disk lifetime on the stellar mass could have strong 
implications for planet formation and migration theories (see e.g. 
[Raymond et al.||2013| |Helled et al.||201 3] [Bandeau et al.|[20l3] 

for recent reviews on this topic), and ultimately result in differ¬ 
ent planetary populations around low-mass and high-mass stars. 
Previous studies have already found hints of protoplanetary disks 


star-forming regions and associations; R14 used this large sam- 

Dahm & Hillenbrand 

2007; Kennedy & Kenyon 

(2009; Roc- 

pie with very low contamination and high completeness levels 

catagliata et al.||2011; 

"ang et al. [2012[ Yasui et a 

. 2014), but 


to derive accurate disk fractions and protoplanetary disk dissipa¬ 
tion timescales. Their results, in good agreement with previous 


estimates (Haisch et al. 2001; Hernandez et al. 

2007b, 2008; Ma- 

majek 2009; 

Fedele et al. 

2010 

Murphy et a 

. 2013) served as 

a benchmark to test the re 

iability and robustness of the sample 


and methods used in R14. 


these studies focused on individual regions and/or were poten¬ 

tially affected by different systematics (e.g., low-number statis¬ 
tics, sensitivity biases, contamination from background giant 
stars, different completeness levels). As a result, a direct com¬ 
parison of their results is difficult and a statistically robust confir¬ 
mation of a possible dependence of the disk lifetime with stellar 
mass has been elusive so far. 
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In this second paper of the series, we study the influence of 
stellar mass on disk evolution using the large and homogeneous 
sample compiled in R14. Section [2] updates the study on R14. 
In Sect. [3] we describe the methodology and criteria used to an¬ 
alyze the disk evolution dependence with stellar mass, and the 
obtained results are shown in Sect. [4] The implications of this 
study are discussed in Sect. [5] 


2. Updates on paper I 

2.1. Description of the sample and bias corrections 

We first make a short summary of the sample compilation and 
completeness corrections in R14. The dataset was obtained by 
combining several stellar population studies of different young 
regions and associations (see Table 1 in R14 for a complete list 
of references). These works used a variety of methods to con¬ 
firm membership depending on stellar masses (e.g. presence of 
lithium in the spectra, accretion, or X-ray properties). Every ob¬ 
ject also has a spectroscopic measurement of its spectral type, 
minimizing the contamination from background giant stars and 
extragalactic sources. Moreover, these studies were specifically 
designed to make a complete census of the stellar content over 
a given mass range and should not be biased towards excess¬ 
bearing stars, for example. We therefore expect the compiled 
sample to be at least representative of (if not complete) the stel¬ 
lar population of each of the considered regions. The covered 
stellar mass spectrum ranges from O-type to late M-type stars 
in all regions. We note that R14 did not re-estimate or discuss 
membership for any target. 

Our main aim in R14 was to derive accurate disk fractions 
for the considered regions via infrared (IR) excesses. This de¬ 
pendence on the availability of IR data could bias the sample, 
and two additional steps were considered to prevent this: 


- spatial correction, i. e., only sources within the field of view 
of IR surveys (Spitzer or WISE, see below) were considered. 
Given that disks are identified via IR measurements, sources 
lacking these data could be incorrectly classified as “disk¬ 
less” (since no IR excess is found), even though there is no 
information to classify them (we do not know whether these 
sources have IR excess or not). We gathered IR data from 
the all-sky WISE survey ( [Wright et al. 2010] ) to study as¬ 
sociations from the SACY sample (Torres et al. |2008| >, and 
therefore they were not affected by this problem. For the rest 
of the regions, however, we used data from different Spitzer 
Space Telescope programs and surveys, which were limited 
in their spatial coverage. In these cases, we checked the posi¬ 
tion of each source and did not considered those outside the 
respective Spitzer map; 

- sensitivity correction, i. e. only sources with predicted pho- 
tospheric level above the sensitivity limit of the correspond¬ 
ing observations were considered. Sources with disks have 
IR excesses, and therefore are easier to detect than naked 
stars (a process known as the Malmquist bias). We therefore 
considered only objects with predicted stellar emission (from 
photospheric models scaled to match the observed fluxes) 
above the corresponding sensitivity limit and wavelength. 
We also took into account the dependence of this limit with 
the corresponding photometric band and survey. For near-IR 
observations, most of the stars fulfill this requirement, and 
therefore the sample is expected to be complete down to mid 
or late M-type stars. In the mid-IR, the sensitivity of the 
observations is more restrictive and the completeness level 
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drops down between late K- and mid M-type stars, depend¬ 
ing on the distance of the region. 

The result of this process is a representative sample of young 
(1-100 Myr) nearby (<500pc) star-forming regions and associ¬ 
ations with minimum biases. 

2.2. Updates in paper I 

[Ribas et ak] ( |2014[ ) identified IR excesses via excess significance, 
defined as: 


f<)hsL'rvccl,.l Fphotosphere,.! 

Xa = - 

& Observed, A 

where F observe(U is the observed flux, F^o^phe,^ is the corre¬ 
sponding model photospheric flux, and crobserved .,1 is the error at 
the corresponding wavelength (A). In R14, bands with ^ > 5 
were considered to harbor excesses. This procedure allowed 
us to compute disk fractions for each region at different wave¬ 
lengths after correcting for sensitivity and spatial completeness 
(see R14 for a description of the method used). 

We updated the study in R14 using new results from the lit¬ 
erature as well as improved methods and calculations. In partic¬ 
ular: 


the association eCha was mistakenly named i] Cha in R14; 
the distance to Serpens has been updated to 415+15 pc from 


Very Long Baseline Array observations of one star (Dzib 
let al.|2010| ); 

we have improved our estimate of the sensitivity complete¬ 
ness limits for each object by taking into account the extinc¬ 
tion. The sensitivity limit should be compared against the 
reddened photospheric value. Reddened photospheres are 
flatter (once normalized) than derredened ones, and hence 
impose looser conditions for the detection limit. As a result, 
the total number of sources within the completeness lim¬ 
its increased significantly, allowing for better statistics and 
slightly varying the disk fractions in R14 (see Appendix); 
the disk identification method in R14 is based on excess sig¬ 
nificances, which depend on the photometric errors (both ob¬ 
servational and systematic). Hence, the disk fractions de¬ 
rived in our study depend strongly on the sensitivity of the 
Spitzer and WISE observations. Smaller errors would result 
in a higher disk detection rate. The photometric data in R14 
were compiled from different programs that used different 
techniques to estimate the photometric errors. Homogeniz¬ 
ing the errors for the whole sample is therefore essential to 
guarantee a constant and coherent disk detection efficiency in 
our heterogeneous dataset. Figure [1] shows the photometric 
uncertainties as a function of apparent luminosity in IRAC3 
and 4 and MIPS1 bands for the sample of 13 star-forming 
regions with Spitzer data. All regions but two (IC348 and 
Ophiuchus) have uncertainties at or below 7% in all IRAC 
3 and 4 bands, and at or below 12% in MIPS1. We inter¬ 
pret the larger uncertainties for IC348 and Ophiuchus as the 
result of the conspicuous and highly variable mid-IR back¬ 
ground observed in both regions at these wavelengths. To 
ensure a homogeneous disk detection rate across the entire 
sample, we rejected Ophiuchus and IC348 for the rest of our 
analysis, and set a conservative constant photometric error of 
7% for IRAC and WISE1, 2, 3 bands, and 12% for MIPS1 
and WISE4 in the x.i calculations. The final sample is com¬ 
prised of 1 809 sources (see Table|A. f[. 
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After these updates, we repeated the analysis in R14 and de¬ 
rived disk fractions in each association. We then fit these disk 
fractions as a function of time using an exponential law of the 
form: 

Ae~' /T + C 


where t is the age of the region (in Myr), and A, r and C are left 
as free parameters. In this simple parametric model, A represents 
the initial fraction of sources with IR excess, r can be interpreted 
as the characteristic timescale of IR excess decay, and C is a 
possible constant level. 

The updated values are presented in the Appendix (Ta¬ 
bles A.l and |A.2| and Figs. |A.l] and |A ,2|. These resu lts are in 
agreement with previou s studies (e.g.|Lada et al.|2000| |Hernan-| 

|dez et alT||2007b[ |2008| |Murphy et al‘||2013[ > and supersede our 

results in R14. 

We note that we have combined data for regions at different 
distances, and therefore the sensitivity limit in spectral type is 
different for each of them. The overall completeness (all regions 
combined) for results with near-IR data (IRAC-WISE1,2,3) is 
around 0.1 M 0 . For mid-IR (MIPS1-WISE4), the sensitivity and 
the completeness limit increases to 0.4 M Q . 


2.3. Impact of new membership estimates and region ages 


New membership estimates are available in the literature for 
some of the regions and associations in this study (e.g. e C’ha, 
TW Hya |Murphy et al.||20l3) |Ducourant et al.||2014| >, which 
could have an impact on our results. However, these works have 
modified the overall members of these associations by a small 
factor and moreover, our results will be still valid as long as the 
considered sample is representative of their stellar population. In 
addition, the original SACY sample has a homogeneous sensitiv¬ 
ity for each region. By including other surveys with completely 
different selection criteria we might add systematics and biases 
difficult to detect and account for. Hence, for consistency with 
R14, we do not update membership or include new objects in 
this work. 

Furthermore, the age estimates for the regions in our study 
are crucial for our analysis. Absolute ages of young associations 
are hard to determine, and some recent works have shown that 
the adopted values may be underestimated by a factor of two 
( |Bell et al.|2013| l. However, the relative age sequence is easier 
to derive and is expected to be more precise (Soderblom et al. 
2013| l. This is actually reinforced by the disk fractions derived 
for these regions: these measurements are independent of the age 
of the regions, yet they clearly decrease once sorted by age (see 
Appendix). We therefore assume the same age values as in R14, 
and we note that our estimates of disk lifetimes may be subject to 
modifications by a scaling factor, depending on the real absolute 
ages. 


3. Stellar masses and ages 

In this section, we study a possible dependence of disk evolu¬ 
tion - as probed by disk fractions - as a function of stellar mass. 
To optimize the disk detection rate, we chose to limit the study 
to regions observed with the Spitzer Space Telescope. Its sen¬ 
sitivity and spatial resolution are significantly better than those 
of WISE, and allow a search for fainter excesses over the photo- 
spheric emission. A total of 1 428 sources from 11 associations 
were observed with Spitzer (see Table[T]i. We note that this value 


Fig. 1 . Ratio of the total uncertainty (observational and systematic) 
over the observed flux as a function of magnitude. IRAC3, IRAC4, 
and MIPS1 bands are shown. Sources from discarded regions (IC348 
and Ophiuchus) are represented as red crosses, the rest of the sources 
as black dots. The black solid line corresponds to the homogenized 
uncertainty levels of 7% (all IRAC bands) and 12% (MIPS1). The 
black dashed line shows the sensitivity limit for each band, as indicated 
by the observing programs. Objects falling beyond this limit in the 
MIPS1 plot are caused by noise (uncertainties in T eff , A v ,...) in the 
corresponding photosphere fit, and were kept in the sample. 



Magnitude 


is smaller than the initial sample, since we are not including as¬ 
sociations studied with WISE data in R14 (i.e., the SACY sam¬ 
ple). The better sensitivity and spatial resolution of the Spitzer 
observations allow us to lower the threshold for excess detection. 
In the following we consider as excess any photometric measure¬ 
ments displaying^- > 3 (instead of;^ > 5 used in Sect. [2]). 

Stellar luminosities were obtained by scaling the best-fit pho- 
tospheric model (BT-Settl models, |Allard et al.||2012) l to the 
dereddened photometric data (see R14) and integrating it over 
the entire stellar spectrum. In this way, the calculation does not 
depend on uncertain bolometric corrections. These luminosities 
were transformed into masses using the mass-luminosity rela¬ 
tionship given by the pre-main sequence isochrones of the cor¬ 
responding age in Bressan et al. ( |2012| ). 

To study the influence of stellar mass in the evolution of pro¬ 
toplanetary disks, we divided the sample into two age bins and 
two mass bins: 

- The age cut separating young and old objects was set to 
3 Myr. This is similar to the typical timescale of disk de¬ 
cay found in Sect|2](see Appendix) for protoplanetary disks. 
A total of three regions (25 Ori, 4 Ori, and Upper Sco) make 
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Table 1. Young nearby star-forming regions included in the study of protoplanetary disk dependence on stellar mass. 


Name 

Age 

(Myr) 

Distance 

(pc) 

Membership and SpT 

Spitzer photometry 

Number of sources 

25 Orionis 

7-10 

330 

(1) 

(1) 

46 

Cha I 

2 

160-165 

(2) (3) (4) 

(3) (4) 

212 

Cha II 

2+2 

178+18 

(5) 

(31) 

47 

CrA 

1-3 

138+16 

(6) (7) (8) (9) 

(32) 

35 

4 Orionis 

4 

400+40 

GO) 

(33) 

114 

Lupus 

1-1.5 

140 - 200 

(11) (12) (13) (14) 

(31) 

217 

NGC 1333 

1 

235+18 

(15) 

(34) 

74 

cr Orionis 

2-3 

440+30 

(16) (17) (18) (19) (20) (21) (22) 

(35) 

104 

Serpens 

2 

415+15 

(15)(23) 

(31) 

142 

Taums 

1-2 

140 

(24) (25) (26) (27) 

(26) (27) 

265 

Upper Sco 

11+2 

140 

(28) (29) (30) 

(29) 

405 


References. (1) Hernandez et al. ( 2007a); (2) Luhman 

2007 ); (3)|Luhman et al. ( 2008a); (4) Luhman & Muench|(|2008); (5) Spezzi et al. (2008 ); 

(6) Neuhauser et al. (20001; (7) Nisini et al. (20051; (8) 

Sicilia-Aguilar et al. (20081; (9) Sicilia-Aguilar et al. (2011); (10) Bayo et al. (2011); (11) 

Krautter et al. (1997); (12) Allen et al. (2007); (13) Comeron (2008); (14) 

Mortier et al. (2011); (15) Winston et al. (2009); (16) Zapatero Osorio 

et al. 

(2002 

; (17) Muzerolle et al. (2003); (18) Barrado y Navascues et al. ( 

2003); (19) Franciosini et al. (2006); (20) Caballero (2007 

i; (21) 

Sacco 

et al. 

(2008 

; (22) Rigliaco et al. (2012); (23) Oliveira et al. (2009); (24) Luhman 

(20041); (25) 

Monin et al. 

(20101); (26) 

Rebul et al. 

(2010 

5 (27) 

Rebull et al. 

(201 1 ); (28) Preibisch et al. (2002); (29) Carpenter et al. 

2006); (30) 

Lodieu et al. 

(2011); (31) 

Evans et al. 

2009 );(32) Peterson et al. 

(201117(33) 

Barrado y Navascues et al. (2007); (34) Gutermuth et al. 

(20081); (35) Luhman et al. (2008b) 


up the “old” sample adding up to 489 sources, whereas the 
rest of the associations (939 sources) make up the “young” 
samplw. 

The mass cut between low- and high-mass objects was set 
to 2 Mq. Although somewhat arbitrary, this threshold is usu¬ 
ally considered as the separation between T Tauri and Herbig 
stars. Given their different physical properties (e.g. radiation 
fields, stellar winds or accretion rates, see |Calvet et al.|2005 ; 


Garcia Lopez et al. |2006[ Hillenbrand|2008| , it is likely that 

disks evolve differently around them. In addition, such a 
mass cut has the advantage of separating the sample so that 
both bins have enough sources for meaningful statistics. 


Therefore, for the rest of the study we classified sources as: 

- “young” if 1 Myr < age < 3 Myr, “old” otherwise 

- “low-mass” if M* < 2 M Q , “high-mass” otherwise 


4. Disk fractions vs stellar mass and time 


The presence of dust around stars is usually detected as IR ex¬ 
cess over the photospheric level: the dust is heated by the stellar 
radiation, and re-radiates with typical temperatures of tens or 
hundreds of Kelvins. However, the amount and shape of this ex¬ 
cess depends strongly on the disk characteristics. Protoplanetary 
disks display excess at short-IR wavelengths (< 10 pm, Willi ams] 
|& Cieza||2011[ ): these are disks that extend from the close re¬ 
gions (< I AU) of the stellar system to hundreds of AUs. How¬ 
ever, the evolution of disks is thought to remove dust from the 
inner regions (e.g. Koepferl et al.||20l3] Alexander et al.||2013 


R14), leaving only mid-IR (24 gm) or even far-IR (> 50 gm) ex¬ 
cesses as in the case of transitional and debris disks ( jWyatt|2008] 
| Williams & Cieza|201 lj ). 


4. 1. Protoplanetary disks 

In this section, we limit the study to short-IR excesses observed 
with IRAC, which offers the best sensitivity to the coolest objects 
of our sample, and therefore allows us to estimate disk fractions 
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over a larger mass range. Unfortunately, IRAC photometry alone 
does not allow a fine classification of disks. We therefore define 
two simple groups: 

- protoplanetary disks are defined here as those sources hav¬ 
ing excesses at wavelengths shorter than 10 pm (either 5.8 
or 8 /.mi, i.e., IRAC3 and/or IRAC4 bands). This definition 
encompasses the classic protoplanetary disks. 

- processed disks are defined as sources without excess ei¬ 
ther at IRAC3 or IRAC4 bands. They could be transitional 
disks, debris disks, or diskless stars, but are indistinguishable 
based on IRAC photometry alone. In all cases their disks (if 
present) have reached a more advanced level of processing 
than the protoplanetary disks defined above. 

We applied the sensitivity completeness correction as de¬ 
fined in Sect. 2 to ensure that the sample is not biased towards 
disk-bearing sources. A total of 1 033 objects met this criteria for 
IRAC4 (the strictest sensitivity limit for IRAC bands). However, 
as pointed in Sect. [2] the considered regions are at different dis¬ 
tances and so the sensitivity limit translates into a different com¬ 
pleteness limit in spectral types for each of them. We computed 
the limiting spectral type for each region, i.e. the latest spectral 
type detected with the corresponding sensitivity limit (see Ta¬ 
ble.^. Given its distance and age, 25 Orionis has a significantly 
earlier limiting spectral type than the rest of the regions, and we 
do not consider it in the analysis. Therefore, we cut the sam¬ 
ple to the most restrictive limiting spectral type (M7, Cha II and 
A Orionis). After this procedure, the resulting IRAC sample is 
comprised of 962 sources. 

We then divided this sample using the mass and age cuts 
described in Sect. [3] The result of this process is shown in the 
form of pie charts in Fig. [2] and in Table [3] Errors for the disk 
fractions were estimated using the bootstraping method further 
described in Sect. |4.3| and are typically less than 6 %. 

After inspection of Fig.[2]some conclusions can be drawn: 

- Protoplanetary disks are significantly more frequent around 
low-mass stars than around high-mass stars, independent of 
the age. 
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Fig. 2. Evolution of protoplanetary disks (red) and processed (yellow) disks for each mass and age bin. The number of sources in each bin is 
shown within brackets. Errors are less than 6 % (see Table [3}. 

Protoplanetary disks evolution 

M„ <2M 0 2M 0 <M, 




- As expected, protoplanetary disk fractions decrease signifi¬ 
cantly at older ages for both low- and high-mass stars. 

4.2. Evolved disks 

As mentioned in Sect. |4.1| IRAC photometry alone does not al¬ 
low a fine classification of processed disks. Photometric mea¬ 
surements at longer wavelengths are required to assess whether 
a source has a transitional or debris disk or if it is a diskless 
star. In an attempt to distinguish these categories, we repeated 
the analysis presented in Sect. |4.1| including MIPS 1 data. The 
24 fi m photometry allows he processed disks to be divided into 
two more precise categories: 

- evolved disks, having no excess at IRAC3 or IRAC4, but 
with an excess at MIPS 1. They include transitional and hot 
debris disks, as well as a small fraction of edge-on disks (e.g 
Huelamo et al. 2010) and circumbinary disks (e.g |Ireland &| 
Kraus|2008| l. 

- diskless stars, displaying no excess at IRAC and MIPS1 
bands. They include genuine diskless stars as well as colder 
(more evolved) debris disks that emit only at far-IR wave¬ 
lengths. 

The same special care was taken to avoid any sensitivity bias, 
and we kept only the sources with predicted photospheric fluxes 
above the sensitivity limit of the observations. MIPS1 observa¬ 
tions are usually less sensitive than IRAC ones. We applied the 
same procedure as in the previous case and excluded the regions 
with significantly worse limiting spectral type (cr Orionis). The 
final sample was limited to spectral types earlier than M4 (the 
corresponding limiting spectral type of CrA). As a result, the 


Protoplanetary (1) 



total number of sources meeting this criteria dropped to 389. 
While smaller and encompassing a narrower mass range, this 
sample allows a more detailed look at the evolution of the inner 
disk via the study of transitional (evolved) disks. 

The sample was then split into the mass and age bins defined 
previously. The results are shown in Fig.[3]and Table[3] We note 
in particular that 

- the results are consistent with those obtained in Sec- 
tion|4.1|for protoplanetary disks within the uncertainties (see 
Sect7[4~3]>; 

- evolved disks are more frequent around high-mass stars than 
around low-mass stars, independently of the age; 

- a roughly constant fraction of evolved disks is found at 
young and old ages. A small difference seems to arise for 
high-mass stars (23 % to 31 %), but it is not statistically sig¬ 
nificant. 


4.3. Significance of the results 

The large number of sources included in the various samples 
in Sects. |4~T| and |4.2| does not guarantee the significance of the 
differences observed in Figs. [2] and [3] and Table [3] Calculating 
uncertainties using standard error propagation or Poisson statis¬ 
tics is neither practical nor trivial in the case of our analysis. The 
computed fractions depend on several different parameters (e.g., 
photometric errors, T eff value and Ay fit, photospheric models 
and adopted isochrones,...) and the age and mass thresholds se¬ 
lected to split the sample could have an impact on the results and 
their interpretation. 
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Fig. 3. Frequencies for protoplanetary (red), evolved (blue) and diskless (gray) objects. The number of sources in each bin is shown in brackets. 
Errors are less than 8 % (see Table[3]l. 

Protoplanetary and evolved disks evolution 

M* < 2 M e 2 M 0 < M t 




Protoplanetary (1) 


75% 




Diskless (33) 



Table 2. Sensitivity limits in spectral type for the IRAC and MIPS1 
samples for each association considered in this study. Ellipses indicate 
that there are no objects in the sample with IRAC/MIPS1 above the cor¬ 
responding sensitivity limit in that region. The adopted limiting spectral 
type for each case is marked with an asterisk. Regions not included in 
the corresponding sample are specified. 


Region 

IRAC sample 

MIPS1 sample 

25 Orionis 

M5 (excluded) 


Cha I 

M9 

M6 

Chall 

M7* 

M5 

CrA 

M8 

M4* 

A Orionis 

M7* 


Lupus 

M9 

M9 

NGC1333 

M8 

M5 

cr Orionis 

M8 

M0 (excluded) 

Serpens 

M9 

M8 

Taurus 

M9 

M6 

Upper Sco 

M9 

M5 


For these reasons, we tentatively estimate uncertainties on 
the disk fractions by performing bootstrapping (1000 iterations) 
randomly varying the T clll photometric fluxes, the young/old 
boundary, the high mass/low mass boundary, and the object 
distances within normal distributions with cr of 50 K, crfi ux , 
0.5 Myr, 0.25 M Q , and error in the region’s distance (or 20 pc 
when no error was available). Uncertainties on the age of the 
sources were not included in the bootstrapping; although abso¬ 
lute ages are hard to determine, relative ages are expected to be 
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more accurate (Soderblom et al. |2013| >. Additionally, our analy¬ 
sis considers only two clearly separated age bins. 

The bootstrapping yielded probability density functions 
(PDFs) of disk fraction for each sample, and should account for 
most uncertainties in the calculations. The standard deviation 
of the PDFs (which are close to normal) provides a reasonable 
estimate of the final uncertainties on the disk fractions and are 
reported in Table [3] 

Overall, the estimated uncertainties suggest that the differ¬ 
ences seen in the disk fractions are statistically significant, and 
that the distributions for high-mass and low-mass stars are dif¬ 
ferent. 

We also checked our results to vary little (a few percentage 
points) when modifying the mass cut between 1.5 and 2.5 M Q , 
and hence our conclusions are robust against the selected mass 
threshold. Given the lack of regions in the 4-8 Myr regime, it is 
not possible to test the impact of a different age cut in detail. Our 
aim is to separate sources which are clearly younger or older, 
and this value is likely not to have a very strong impact on our 
study. For this reason, our results are weakly dependent on age 
uncertainties, as long as relative ages are properly estimated (see 
discussion in Sect. |2.3| l. 

4.4. Disk fractions and stellar temperature (mass) 

Our definition of protoplanetary, evolved and diskless sources 
depends on the wavelength at which excess over the photosphere 
is found. Because the temperature distribution in the circumstel- 
lar disk depends directly on the stellar effective temperature and 
luminosity, the classification probably suffers from a dependence 
on stellar temperature. Infrared emission at a certain wavelength 
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Table 3. Results from the analysis of disk frequencies as a function of age and mass. See also Figs. [ 2 ] and [5] Errors are one standard deviation 
from the computed PDFs. 


Disk type 

Young+Low-mass [%] 

Young+High-mass [%] 

Old+Low-mass [%] 

Old+High-mass [%] 

Protoplanetary disks (Section 4.1 1 

Protoplanetary (IRAC4 detection limit) 

63±2 

38±6 

17±2 

9+4 

-2 

Processed (IRAC4 detection limit) 

37±2 

62±6 

83±2 

98/2 

Protoplanetary and evolved disks (Section 4.2\ 

Protoplanetary (MIPS1 detection limit) 

66±4 

36±8 

8±2 

9+4 

-2 

Evolved (MIPS 1 detection limit) 

9±2 

31 ±7 

8±3 

23±6 

Diskless (MIPS1 detection limit) 

25±3 

33±8 

84±3 

75±7 


likely originates from different regions of the disk for hot (mas¬ 
sive) and cool (low-mass) stars, and may even arise from differ¬ 
ent disk types. 

This effect is probably not critical for our analysis of pro¬ 
toplanetary sources (Sect. |4T| . These dense and massive disks 
have complex physical structures, and the IR excess at a cer¬ 
tain wavelength is emitted from relatively large areas of the disk. 
Moreover, the protoplanetary definition is based on the presence 
of excess over a relatively wide wavelength range, ensuring that 
the corresponding emission comes from broad areas of the disk 
regardless of the stellar mass. 

On the other hand, our definition of evolved disks relies pri¬ 
marily on IR excess at 24gm (MIPS1). The area of the disk 
probed at 24 pm is expected to vary significantly between mas¬ 
sive and low-mass stars, and might make the comparison of the 
excess fractions and their interpretation more difficult. 

To test how much this dependence might affect our results, 
we computed a grid of radiative transfer models of circumstellar 
disks using the code MCFOST ( |Pinte et al.|2006| l. We found that 
synthetic SEDs classified as evolved according to our criterion 
always corresponded to higher levels of disk processing regard¬ 
less of the star’s temperature. These simple simulations suggest 
that our classification scheme is efficient and robust enough to al¬ 
low a comparison of the results obtained for high and low-mass 
stars. 

Finally, a comparison of the evolution of evolved disks can 
be performed within the same mass bin, removing most of the 
previous dependency on stellar temperature. Figure [4] shows 
the PDFs of the "evolved" disk fractions obtained in Sect. IQ 
for high-mass stars on the one hand and low-mass stars on the 
other hand. The PDFs are completely compatible both for low- 
mass and high-mass stars, pointing to a roughly constant level 
of evolved disks within the uncertainties at any age regardless of 
stellar mass. 

5. Impact of stellar mass on disk evolution 

5.1. Dependence of protoplanetary disk lifetimes on stellar 
mass. 

Previous observational studies had already found lower proto¬ 
planetary disk fractions around high-mass stars and suggested 
that their protoplanetary disks disperse earlier and/or evolve 
faster (see e.g Williams & Cieza [2011 for a review ). A number 
of these studies focused on individual regions, making it hard 
to confirm a global dependence of protoplanetary disk evolution 
on stellar mass independently of the environment and method 


used to identify the disks (e.g 

^ada et aL|2006, 

Dahm & Hillen- 

brand 2007 1 Currie et al.|2009 

Carpenter et al. 

2006 

in IC348, 

NGC 2362 and Upper Sco). 

Cennedy & Kenyon 

2009) per- 


Fig. 4. Comparison of the obtained PDFs for evolved sources as a 
function of age and for the two stellar mass bins. 




formed an original and coherent study of disk frequencies in nine 
young associations by combining optical spectra and IR photom¬ 
etry. In contrast to our study, they found only a weak depen¬ 
dence of disk dispersal on stellar mass. The ages used in their 
study have been improved and sometimes changed significantly. 
For example, the age of Upper Sco was recently re-estimated to 
11 Myr|Pecaut et aL] ( |2012j l from an original 5 Myr ( |Preibisch~&| 
Zinnecker|1999| l, resulting in a noisier distribution of disk frac¬ 
tion as a function of age in their study. 

Our results demonstrate quantitatively that protoplanetary 
disks dissipate significantly faster and/or earlier around high- 
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Fig. 5. Evolution of protoplanetary disks as a function of mass, as 
derived in this study. 
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mass stars, as we find between 10 and 30 % more protoplanetary 
disks around low-mass stars regardless of the age. 

Our data additionally show that the disk fraction for high- 
mass stars falls to ~()% anytime between 3 and 11 Myr, and 
suggest that protoplanetary disks disperse up to two times faster 
around high-mass stars than around low-mass ones. More data 
for regions with intermediate ages are needed to better constrain 
this ratio. 

Interestingly, other studies found a higher (80-85 %) proto¬ 
planetary disk fraction and no dependence with stellar mass in 
the younger Trapezium cluster (< 1 Myr, [ Lada et al . 12000) and 
NGC 2024 cluster (<0.4-1.4xl0 5 yr, |Hais(Feral.|2000| )~Com- 


bined with the present study, these results suggest that 


- the disk fractions before 1 Myr are in the continuity of the 
results presented here; 

- the dependence of disk fractions on stellar mass appears af¬ 
ter the first Myr of the star’s lifetime, and persists until no 
protoplanetary disks are found around ~ 10 Myr 

Based on these two suggestions, we propose a scheme of disk 
evolution as a function of stellar mass (see Fig. [5]). The timescale 
of protoplanetary disk formation and evolution is mostly inde¬ 
pendent of stellar mass until the first Myr. At 3 Myr, significant 
differences are already observable between the evolution of pro¬ 
toplanetary disks around low- and high-mass stars. By ~ 10 Myr, 
no protoplanetary disks are found around high-mass stars, but 
*10-15% of low-mass stars still harbor a disk. Stronger radi- 


ation fields and higher accretion rates of high-mass stars ( 

Cal- 

vet et al. 2005; Garcia Lopez et al. 2006 1 Hillenbrand 2008 

are 


likely to affect the evolution of protoplanetary disks even during 
the first Myrs, fastening the processes responsible of disk disper¬ 
sal ( [Alexander et al.|2013| l. 


5.2. A constant level of evolved disks 

We find the fraction of evolved sources to remain constant with 
age both for low- and high-mass stars: 5-15 % of low-mass stars 
display this kind of excess, and the fraction increases to 20- 
30% for high-mass objects (see Fig. 4]). Again, the dependence 
with stellar mass of the disk radii probed at a certain wavelength 
should be consider and it could, in fact, account for the higher 
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fraction of evolved disks around high-mass stars. Nevertheless, 
the obtained values at these ages can be established to be within 


5-30%. 


The evolved disk definition proposed here includes transi¬ 
tional disks, hot debris disks, edge-on disks and circumbinary 
disks. Unfortunately, the classification does not allow us to dis¬ 
criminate between these four types of disks. We nevertheless 
note that simple geometric considerations suggest that edge-on 
disks must be relatively rare contribute only moderately to the 
population of evolved disks presented here. Hot debris disks are 
expected to appear mostly at later stages of disk evolution and 
must populate mostly the older bin of our analysis. Circumbi¬ 
nary disks are expected to be common given the multiplicity fre¬ 
quency among stars (e.g|Bouy||2( 


y 2011) , but only the closest sys¬ 
tems would mimic the SED of evolved disks defined here (e.g as 
in |Ire land & Kraus|2008 Nguyen etal.|2012|l. The evolved disks 


in the young (< 3 Myr) bin are therefore most likely primarily 
made of transitional and circumbinary disks, and our analysis 
suggests that these two types of disks can be found even during 
the very first Myrs of stellar life. 


Finally, we also note that the higher fraction of evolved 
disks - which include circumbinary disks - around high-mass 
stars is also consistent with their higher multiplicity frequency 
( [Duchene & Kraus|2013) l. 


5.3. Transitional disks 


We stress that the current data in our sample cannot be used to 
isolate transitional disks among evolved disks. With this limita¬ 
tion in mind, it is nevertheless interesting to compare the prop¬ 
erties of evolved disks of the present study to the properties of 
transitional disks as reported in the literature. 

Transitional disks are a particularly interesting kind of cir- 
cumstellar disks that harbor a gap and/or cavity ( [Andrews et al.) 
201 1[ >. Various mechanisms have been proposed to explain 
the existence of these gaps and cavities. They can be pro¬ 
duced by the formation of planets ( |Espaillat et al.||20T4| ), grain 
growth, photoevaporation or the formation of low-mass compan- 
ions (|Birnstiel et al.||2012| [Alexander et al.|[2013| |Papaloizou| 
|et al.|2007| ). The connection between transitional disks and plan¬ 
ets - as demonstrated by several resolved images of planetary 
mass companions within the gap or cavity (e.g [H ue lamp et al. 
2011 |Kraus & Ireland||2012| |Quanz et al.||2013[ |Close et al. 
2014) - makes them a cornerstone for theories of planet forma¬ 
tion and evolution. 


First, we note that the fraction of evolved disks derived by 
our analysis (between 5 and 30%) is consistent with transitional 


disk fractions reported in the literature (e.g. Skrutskie et al. 

o 

On 

On 

Kenyon & Hartmann 1995 Andrews & Williams 2005; 

Cieza 

et al. 2007. Men'll et aL|2010). 


Second, previous studies have established the duration of the 
transition phase to be s: 1 Myr (Alexander et al. 2013), which 


is in agreement with the combined 5-30% of evolved sources 
and the « 10 Myr disk lifetime, if one assumes that all disks go 
through a transitional phase. Within the planet formation sce¬ 
nario, only Jovian planets are likely to produce a gap in their 


protoplanetary disk (pin & Papaloizou|1986 

Marsh & Mahoney 

1992; [Nelson et al. 2000. Calvet et al. 2002 

Rice et al. 2003), 


and so it is also plausible that not all planetary systems go 
through a transitional phase during their formation. If that is the 
case, then the estimated ~ 1 Myr duration of this phase should be 
regarded as a lower limit. 
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5.4. Implications for planet formation and giant planet 
migration 

To date, there are two main theories explaining the formation of 
planets: the core accretion model (CA) and the disk instability 
model (DI) (see |Helled et al.||2013) [Raymond et al.||2013| for 
recent reviews on this topic): 

- In the CA scenario, a solid planetary embryo is formed by 
continuum accretion of planetesimals onto this embryo. Ter¬ 
restrial planets are likely formed in this way. If the solid core 
becomes massive enough, it will also start accreting gas at 
an increasing rate until no more gas can be obtained from 
the disk (either because the disk has dissipated or because 
the planet has opened a gap in it). The result of this process 
is a gas giant planet. 

- In the DI model, a very massive disk is destabilized by its 
own gravity, leading to the collapse of a region of the disk 
and potential formation of a giant planet or brown dwarf. 

So far, the CA model seems to better reproduce several ob¬ 
served properties of the known exoplanets. For example, cor¬ 
relations are known to exist, both between the stellar metallic- 
ity and giant planet occurrence, and between stellar and plane¬ 
tary metallicities (e.g. lohnson et al. 2010j ). These correlations 
can be explained by the CA model, in which higher metallic¬ 
ities result in more efficient planet formation. However, some 
known exoplanetary systems cannot be explained within the CA 
frame (e.g. the system around HR 8799, |Marois et al.||2Q08] >, 
suggesting that DI may be responsible for the formation of some 
gas giants at very large radii from their host stars and/or around 
low-metallicity stars. It could also explain the existence of some 
brown dwarfs orbiting other stars. The global planet formation 
process may be a combination of these two scenarios. 

If CA is the dominant mechanism, then planets are likely 
to form below the ice-line ( |Ida & Lin|2004a| >: water can freeze 
into ice around the grains of dust, and increase the efficiency of 
their aggregation to form the required planetary embryo. On the 
other hand, formation at very large radii from the host star is 
also unlikely, since the solid accretion rate of planetary embryos 
decreases with the distance to the star. As an example, the op¬ 
timal zone to form planets for a solar-type star is estimated to 
be between 5-10 AU ( jHelled et al.||2013] >. Surprisingly, a sig¬ 
nificant number of Jovian-like planets are found very close to 
their stars, and are referred to as hot Jupiters. This particular 
subgroup of exoplanets suggests that a fraction of these gas gi¬ 
ants migrate inwards after their formation: once a gap is opened 
by the planet, it gets locked to the disk and exchanges angular 
momentum with it, causing the planet to move closer to its host 
star (Type-II migration, Papaloizou et a l.[2007) l. Other migration 
processes such as planet-planet scattering, Kozai resonances, or 
similar dynamic mechanisms (e.g. |Kozai|[l962[ Wu & Murray | 
2003 [ Naoz et aT7||2011)> are likely responsible for the fraction 


of hot Jupiters whose spin is found to be misaligned with their 
planet-star orbits dTriaud et al. 120101 INarita et al .120 1 01 IMortonl 
|& Johnson|2011| l. 

Within this theoretical framework, the obtained dependence 
of stellar mass and disk lifetime would have two main implica¬ 
tions: 

1. Observational results tentatively suggest that massive stars 
are more likely to host gas giants than their low-mass coun¬ 
terparts (Johnson et al. 2010| >. This trend would be the re¬ 
sult of two competing factors: high-mass stars have more 
massive disks (e.g. Andrews et al. 2013), and so are more 


Fig. 6. Orbital semi-major axes as a function of the mass of the host 
star for the known exoplanet population. Red circles show hot Jupiters 
(M > lMj, period < 10days), black circles represent the rest of the 
sample. The size of the symbols scales with the planet mass. We also 
show the estimated stellar mass radii of main-sequence stars ( jSiess et ah] 
[2000) and the stellar radii of the host stars. A clear gap of hot Jupiters 
is seen for M* > 1.5 M 0 . 

it 
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likely to form gas giants. On the other hand, the present 
study shows that massive stars disperse their disks earlier. 
Under some circumstances, the formation of giant planets 
could even last for some Myr ( Helled et al.||2013j ). There¬ 
fore, a fast dispersing disk may not live long enough to form 
gas giants. If the observed trend (massive stars being more 
likely to host giant planets) is real, then either a strong corre¬ 
lation between the stellar mass and the ability to form giant 
planets exists to compensate for the shorter disk lifetime, or 
giant planet formation is faster than the disk dispersal time 
regardless of the stellar mass, or most of the planet formation 
process happens in the first million years. 

2. Type-II migration is active while the gaseous disk is present 
and the planet can exchange angular momentum with it. The 
distribution of orbital semi-major axes produced by these mi¬ 
gration mechanisms is controlled by the ratio of the disk de¬ 
pletion timescale to the viscous diffusion timescale at the 
distance where Jovian planets form (Armitage et al. 2002| 
|Trilling et al.||200~2[ |Ida & Liii||2004 r|l. A fester dispersal 
of protoplanetary disks around high-mass stars would result 
in a shorter disk depletion timescale, and therefore a shorter 
migration period ( Burkert & Ida|2007 |Currie|2009| l. Our re¬ 
sults provide a robust confirmation of a faster disk dispersal 
around high-mass stars: the disk lifetime could be more than 
two times shorter for massive stars, reducing significantly the 
migration period available for Jovian planets and, as a result, 
a smaller number of hot Jupiters would be present around 
these stars. 


Following Kennedy & Kenyon ( 2009| ), point 2 can be linked 
with the known population of exoplanets. Figure [6] shows the 
orbital semi-major axis distribution as a function of the mass 
of the host star for the known exoplanet^] (planets around pul¬ 
sars are not plotted since they represent a completely different 
type of object). Considering as hot Jupiters the planets with 
masses above 1 My and periods < 10 days, a clear lack of these 


as in the www. exoplanets. eu database, August 6 th 2014 
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objects around massive stars (> 1.5 M Q ) is easily identified. This 
paucity of close Jovian planets around massive stars is still a mat¬ 
ter of intense debate. Tidal interaction between massive planets 


Livio 2009; Villaver et al. [2014 

. Some authors report a real de- 

pendence on stellar mass (e.g. 

Johnson et al. 2011), and have 


even shown hints of it not being the result of planet engulfing 
caused by the post-main-sequence evolution of the star ( jJohn-[ 
son et al.|2007[>. On the other hand, other authors suggested that 


such claims were made on incorrectly assigned stellar masses 
(Lloyd 201 l|l, and that tidal destruction is responsible for the 


2013|). By also plotting the 


observed trend (Schlaufman & Winn 

radii of the planet host stars in Fig. [6 it becomes clear that most 
low-mass stars (< 1.5 — 2 M 0 ) have radii matching those pre¬ 
dicted from main-sequence models (Siess et al. |2000 >. For more 
massive stars, the stellar radii increase owing to the post-main- 
sequence phase, overlapping with the loci of the hot Jupiters 
semi-major axes. It is not possible to tell if the massive stars 
in Fig.[6]have engulfed any close gas giants, but it becomes clear 
that the post-main-sequence phase of these stars could account 
for a fraction of the observed paucity of hot Jupiters around 
them. Assessing the origin of this lack of close gas giants is 
beyond the scope of this study, but the faster disk dissipation 
around high-mass stars could contribute to it as well. 

This discussion may also apply to some brown dwarfs. Al¬ 
though it is thought that most of these objects form via differ¬ 
ent mechanisms (Chabrier et al. |2014) >, the existence of binary 
systems including a brown dwarf companion suggests that some 
of them may form in a circumstellar disk around more massive 
stars. If this is the case, they are very likely formed in similar 
ways to giant planets, and the same conclusions could be applied 
to this particular subgroup of brown dwarfs. 

Finally, the implications of the faster disk dispersal around 
high-mass stars for terrestrial planets are harder to determine. As 
previously mentioned, their formation follows the initial stage of 
the CA model, accreting planetesimals to form a planetary em¬ 
bryo and eventually a rocky planet, which will not become a 
gas giant if the planet cannot accrete enough gas for any reason. 
Flowever, this process is faster than gas giant formation: earth- 
mass planets can be formed in 1 Myr ( |Raymond et al.|2013) >, and 
so the disk dissipation dependence with stellar mass may not be 
as relevant as for giant planet formation. Moreover, terrestrial 
planets do not undergo Type-II migration, and therefore their or¬ 
bits are less dependent on the disk properties. However, their 
orbits can be strongly affected by the migration of giant planets 
in the same system. If the migration history of gas giants has 
a dependence with stellar mass, it could indirectly produce dif¬ 
ferent orbits for terrestrial planets too. The identification of any 
possible trend of this effect requires further modeling and study 
of the known planetary population. 


6. Conclusion 

We have studied the disk lifetime dependence on stellar mass 
using the sample of YSOs compiled by Ribas et al. 1(2014} . Af¬ 
ter updating their results, we have divided the sample into low- 
or high-mass stars (boundary set to 2 M G ) and young or old re¬ 
gions (boundary set to 3 Myr). We then study the fraction of 
protoplanetary and evolved disks as a function of age and mass. 
The large number of sources in the sample allows us to con¬ 
firm the existence of a dependence of disk lifetime on the stellar 
mass: high-mass stars disperse their disks up to two times faster 
and/or earlier than low-mass ones. We also find that the fraction 
of evolved disks (including transitional, circumbinary, and hot 
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debris disks) remains roughly constant (5-30 %) during the first 
10 Myr of the stellar life. The faster dispersal of protoplanetary 
disks around high-mass stars could have important implications 
for giant planet formation and migration, and may contribute to 
the apparent lack of hot Jupiters around these stars. 
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Table A.l. Fraction of sources with excess for each region in the three wavelength regimes (update of Table 4 in R14) 


Name 

Fraction i RA c(%) 

Fraction sho rt(%) 

Fraction intermediate (%) 

Fraction i ons (%) 

25 Orionis 

12 + 7 [26] 

8 + 5 [26] 

12 + 7 [25] 

...[0] 

Chal 

51 +4 [117] 

36 + 5 [109] 

52 + 5 [106] 

71+6 [52] 

Cha II 

75 + 8 [28] 

50 + 9 [28] 

85 + 9 [20] 

82 + 7 [28] 

CrA 

40+ 10 [19] 

11 + 8 [19] 

50+ 10 [16] 

...[6] 

IC 348 

38 + 3 [253] 

17 + 2 [262] 

41 + 3 [223] 

81+8 [27] 

/I-Orionis 

23 + 7 [43] 

10 + 4 [51] 

23 + 6 [43] 

...[0] 

Lupus 

53 + 5 [85] 

27 + 5 [86] 

53 + 5 [85] 

57 + 6 [58] 

NGC 1333 

64 + 5 [73] 

38 + 6 [73] 

67 + 6 [70] 

88+8 [17] 

Ophiuchus 

23 + 3 [248] 

13 + 2 [253] 

23 + 3 [238] 

42 + 5 [113] 

cr-Orionis 

34 + 6 [71] 

14 + 4 [74] 

34 + 6 [71] 

...[0] 

Serpens 

57+ 4 [129] 

35 + 4 [131] 

58 + 5 [128] 

57+7 [56] 

Taurus 

61 +4 [202] 

39 + 3 [214] 

62 + 3 [197] 

76 + 6 [49] 

Upper Sco 

9 + 2 [250] 

4 + 1 [250] 

10 ± 2 [232] 

14 + 3 [186] 

AB Dor 

4 + 2 [77] 

4 + 2 [77] 

1 + 2 [77] 

0 + 1 [50] 

Argus 

0 + 2 [45] 

0 + 1 [44] 

0 + 2 [44] 

...[3] 

P Pic 

8+4 [50] 

4 + 3 [50] 

6 + 3 [49] 

14 + 5 [43] 

Carina 

4 + 3 [28] 

0 + 1 [27] 

4 + 3 [28] 

...[6] 

Columba 

2 + 2 [57] 

0.0 + 0.9 [57] 

2 ± 2 [57] 

10 + 8 [20] 

e Cha 

21+8 [24] 

17 + 7 [24] 

17 ± 8 [24] 

20+ 10 [13] 

Octantis 

0 + 3 [17] 

0 + 2 [17] 

0 + 3 [17] 

...[1] 

Tuc-Hor 

6 + 3 [48] 

6 + 3 [48] 

0 + 1 [48] 

2 + 3 [41] 

TWHya 

20 + 9 [20] 

0 + 3 [20] 

20 + 9 [20] 

20+ 10 [14] 


Notes. The number of sources used is indicated within brackets. Sub-indexes reference fractions at different wavelength ranges, as defined in R14 


Table A.2. Values and errors obtained from fitting and exponential decay to disk fractions with time. 


Wavelength range 

A 

(%/Myr) 

T 

(Myr) 

C 

(%) 

Short 

60+ 10 

2.7 + 0.7 

1.1 +0.9 

Intermediate 

90+ 10 

3.3+ 0.6 

1 + 1 

Long (primordial) 

84 + 6 

4.4 + 0.5 

-0.1 +0.4 

Long 

95 + 8 

5.3+ 0.9 

0 + 1 


Notes. The long (primordial) case corresponds to red circles in Fig.|A.2| 
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Fig. A.l. Upper left panel: comparison of the disk fractions obtained in this study (circles) with those from (Hernandez e t al.| 2007b [20081 
(squares) for regions in common. Upper right and lower panels: disk fractions for all regions in the short, intermediate and long wavelength 
ranges. The best-fit exponential law is over-plotted (red line). 



Age (Myr) 


Age (Myr) 


Fig. A.2. Fraction of primordial (red circles) and evolved (blue squares) disks as a function of age. The best-fit exponential law for the primordial 
disk percentages is also shown (red dashed line). 
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